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Abstract 

Objective: Osteoarthritis (OA) is a multi-factorial disease and genetic factor is one of the important etiologic risk 
factors. Various genetic polymorphisms have been elucidated that they might be associated with OA. Recently, 
several studies have shown an association between Growth Differentiation Factor 5{GDF5) polymorphism and knee 
OA. However, the role of genetic predisposing factor in each ethnic group cannot be replicated to all, with 
conflicting data in the literatures. Therefore, the aim of this study was to investigate the association between GDF5 
polymorphism and knee OA in Thai population. 

Materials and Methods: One hundred and ninety three patients aged 54-88 years who attended Ramathibodi 
Hospital were enrolled. Ninety cases with knee OA according to American College of Rheumatology criteria and 
one hundred and three cases in control group gave informed consent. Blood sample (5 ml) were collected for 
identification of GDF5 {rsl 43383) single nucleotide polymorphism by PCR/RFLP according to a standard protocol. 
This study protocol was approved by the Ethics Committee on human experimentation of Ramathibodi Hospital 
Faculty of Medicine, Mahidol University. Odds ratios (OR) and 95% confidence intervals were calculated for the risk 
of knee OA by genotype (TT, TC and CC) and allele (T/C) analyses. 

Results: The baseline characteristics between two groups including job, smoking and activity were not different, 
except age and BMI. The entire cases and controls were in Hardy-Weinberg equilibrium (p > 0.05). The OA knee group 
(n = 90) had genotypic figure which has shown by TT 42.2% (n = 38), TC 45.6% (n = 41 ) and CC 1 2% (n = 1 1 ), whereas 
the control group (n = 103) revealed TT 32% (n = 33), TC 45.6% (n = 47), and CC 22.3% (n = 23), respectively. Genotypic 
TT increased risk of knee OA as compared to CC [OR = 2.41 (P = 0.04, 95%CI = 1.02-5.67)]. In the allele analysis, the T 
allele was found to be significantly associated with knee OA [OR = 1 .53 (P = 0.043, 95%CI = 1 .01-2.30)]. 

Conclusion: These data suggested that GDF5 polymorphism has an association with knee OA in Thai ethnic. This 
finding also supports the hypothesis that OA has an important genetic component in its etiology, and GDF5 
protein might play important role in the pathophysiology of the disease. 
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Introduction 

It is widely believed that osteoarthritis develops from an 
imbalance between anabolic and catabolic processes or 
homeostasis of cartilage metabolism [1,2]. The etiology 
of this disease is related to genetic association [3]. 
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Recently, several studies have demonstrated the poly- 
morphism in many genes which might be related to the 
pathogenesis of osteoarthritis [4-14]. Growth Differentia- 
tion Factor 5 or GDFS gene regulates the expression of 
the GDFS protein which is closely related to BMP and 
is a member of TGF-beta superfamily [15]. It has a role 
in regulation of the chondrogenesis and the defect in 
this gene might be correlated to the abnormal joint 
development [15]. It has been reported in animal study 
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that GDF5 knockout mice develops knee joint anomaly 
[16]. Moreover, it has been shown that he polymorph- 
ism in GDF5 gene is related with low expression of the 
GDFS protein in knee joint [17]. 

The large scale analysis has shown that the association 
of this polymorphism (rsl43383) in promoter area 
might be a risk factor of the osteoarthritis [18]. In addi- 
tion, studies of this genetic variant in China and Japan 
have shown the association of T allele and knee OA 
[19]. However, there is a report from Greece which had 
an inconsistent result and found no an association [20]. 
The genetic susceptibility of the disease in different eth- 
nic cannot be applied to others because each ethnic has 
a different genetic background. There is a gap of infor- 
mation about this polymorphism in Thai population, 
therefore the objective of this study is to determine the 
association of the SNP rsl43383 in GDFSgene and knee 
OA in Thai population. 

Patients and Methods 

Subjects 

Our study was approved by the Ethics Committee of the 
Faculty of Medicine, Ramathibodi Hospital, Mahidol 
University, Bangkok, Thailand. All patients recruited in 
the study were Thai by nationality and had ancestors 
settled in Thailand for at least three generations. In 
total, 90 patients with knee OA who underwent total 
knee arthroplasty (TKA) and 103 patients without knee 
OA were enrolled from Department of Orthopedics, 
Ramathibodi Hospital. Informed consent was performed 
after the purpose of the research project and it had 
been clearly explained to the patients. The diagnosis of 
knee OA was based on the American College of Rheu- 
matology criteria [21]. Both OA and control groups 
were interviewed to obtain demographic data and all of 
established risk factors. Thereafter, standard weight- 
bearing antero-posterior and lateral view of knee radio- 
graphs were taken to confirm the diagnosis of OA by 
Kellgren and Lawrence scores (KL scores) [22]. 

Laboratory technique 

PCR-RFLP for BsiEl restriction site was used for SNP in GDF- 
5 identification 

Firstly, 5 ml peripheral blood sample was collected from 
patient using ethylenediamine tetraacetic acid as an 
anticoagulant and processed for SNP analysis. Genomic 
DNA was extracted from buffy coat leukocytes using the 
standard phenol-chloroform method. PCR primers to 
amplify the promoter area of GDF5 gene were designed 
by the Primer-3 web-based tool [23]: GATTTTTTCT- 
GAGCACCTGCAGG (forward) and GTGTGTGTT 
TGTATCCAG (reverse). 50 ul PCR mixture contained 
100 ng of genomic DNA, 20 pmol of each primer, 0.2 u 
M of each dNTP, 1 unit of Taq DNA polymerase 



(AmpliTaq®, Applied Biosystem, Foster City, CA), 3.0 
mM MgCl2 in 10 x PCR buffer containing 10 mmol of 
Tri-HCl pH 9.0, 10 mmol KC1 and 0.1% Triton X-100 
(Invitrogen, Carlsbad, CA). PCR reaction was started 
with an initial denaturation at 95°C for 5 min, followed 
by 35 cycles of amplification in a thermocycler (PCR 
Sprint, Thermofisher, Waltham, MA) with denaturation 
at 94°C for 1 min, annealing at 58°C for 1 min, exten- 
sion at 72°C for 1 min, and final extension at 72°C for 
10 min. Then, 10 \iL of PCR product was incubated at 
37°C with 3 units of BsiEI for 4 hours under the manu- 
facturer's recommended conditions (New England Bio- 
labs, Ipswich, MA). The digested product was 
electrophoresed on 2% agarose gel with ethidium bro- 
mide staining before being visualized on a UV transillu- 
minator (Figure 1). The expected fragment length was 
104 and 230 bp in CC, 104, 230, and 344 bp in TC, and 
344 bp in TT genotypes, respectively. 

Statistical analyses 

Analysis of demographic data was performed by Excel 
2007 (Microsoft® Excel®). The unpaired t-test was used 
for continuous data and Chi-square was used for catego- 
rical data. Allele frequencies, Odds ratio and the prob- 
ability for Hardy- Weinberg Equilibrium were estimated 
and analyzed as explained by following website; http:// 
ihg.gsf.de/cgi-bin/hw/hwal .pi. 

Results 

In our study, the average age of patients in the knee OA 
group was significant older than in the control. BMI 
which is one of the risk factor of knee OA was also 
slightly higher in the OA patient group. However, other 



























Figure 1 shows the bands of 
electrophoresis gel: A represents TT 
represents TC. 


each genotype in the 
, B represents CC and C 
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possible risk factors which might be related to OA as 
well as patient life style were comparable between these 
two groups. The baseline characteristics of the patients 
are shown in Table 1. The prevalence of allele T and C 
in our sample was normally distributed according to the 
Hardy- Weinberg Equilibrium (p-value < 0.05). The com- 
parison of the number and odd ratio of genotype TT, 
TC, TT + TC with CC of GDFS polymorphism 
(rsl43383) in the promoter area between case and con- 
trol is demonstrated in Table 2. When analyzing by 
allele (shown in table 3), the odds ratio of T allele was 
1.53 (95%C.I. = 1.01-2.31), which delineated the increase 
in risk to develop knee OA. Our result shown that this 
polymorphism, T/C, is inherited by autosomal recessive 
manner as the TT genotype increases risk of the disease 
significantly (OR = 2.41, 95%C.I. = 1.02-5.67), whereas 
the TC has no significant difference. 

Discussion and Conclusion 

Our study has shown the statistically significant associa- 
tion between polymorphism in promoter of GDFS gene 
(rsl43383) and knee osteoarthritis in Thai population. 
We found that T allele in GDFS polymorphism was a 
significant risk factor for knee osteoarthritis with Odds 
ratio 1.53 (95%C.I. = 1.01-2.31). When analysed by gen- 
otype, it was found that the TT genotype increased risk 
of knee OA, whereas TC has not shown any correlation. 
Regarding these finding, it suggests that this poly- 
morphic gene might be expressed as the autosomal 
recessive type. The GDFS gene is located on the chro- 
mosome 20qll.2 and it regulates the expression of 
GDFS protein. It is categorised in the class of bone mor- 
phogenetic protein (BMP) [15]. It involves in develop- 
ment of bone and cartilage, particularly in endochondral 
ossification process [24,25]. The mutation of GDFS 
associates with generalised osteoarthritis and skeletal- 
related congenital diseases. Acromesomelic chondrodys- 
plasia; Grebe type and Hunter-Thompson type is 



Table 1 Shows the base line characteristics between 
cases and controls 



Variables 


Cases 
(n = 103) 


Control 
(n = 103) 


p-value 


Age (years), mean (SD) 


68.46(10.0) 


59.25(9.0) 


< 0.01* 


Onset of disease (years), mean (SD) 


61.46(9.06) 


N/A 


N/A 


Female (%) 


79(87.25) 


93(90.29) 


0.7 


BMI (kg/m 2 ), mean (SD) 


26.66(4.0) 


24.57(5.0) 


< 0.01* 


History of labor work (%) 


25(27.78) 


30(29.13) 


0.20 


Regular exercise (%) 


10(11.11) 


13(12.62) 


0.15 


Kneeling activity (%) 


37(41.11) 


48(46.60) 


0.07 


Previous any fracture (%) 


9(10.0) 


13(12.62) 


0.15 


Smoking (%) 


2(2.22) 


4(3.88) 


0.10 



Note: ^statistically significant, SD = standard deviation, N/A = not applicable 



Table 2 Shows the association between each genotype 
and risk of knee OA 



Genotype 


Case 
(n = 90) 


Control 
(n = 103) 


Odds ratio 
(95%CI) 


p-value 


TT 


38 


33 


2.41(1.022-5.670) 


0.04* 


TC 


41 


47 


1.82(0.794-4.190) 


0.15 


^ + TC 


79 


80 


2.07(0.944-4.517) 


0.06 


CC 


11 


23 


1 


N/A 



(CC was used for calculation Odd ratio of TT, TC and TT+TC) 



Notes: ^statistically significant, N/A = not applicable 

autosomal recessive form which is rare hereditary skele- 
tal disorders. They are presented by short status, abnor- 
mal limbs development. Brachydactyly type Al, A2 and 
C and symphalangism proximal syndrome are inherited 
malformation presented with abnormal morphology of 
hand and finger. Multiple synostoses syndrome type 2 is 
an autosomal dominant condition characterised by pro- 
gressive joint fusions and progressive conductive deaf- 
ness. Du Pan syndrome is a rare autosomal recessive; 
the patient presents with absence of the fibulae and 
severe acromesomelic limb shortening with small, non- 
functional toes. These abnormal conditions have been 
reported that they have a mutation in GDF5 gene, thus 
they support that GDF5 plays a role in skeletal develop- 
ment [26,27]. Recently, it has been reported that GDF5 
deficiency mice had the delay fracture healing process 
that impaired the cartilaginous matrix deposition in the 
callus and reduced callus cross-sectional area [28]. 

The functional study of this polymorphism has been 
demonstrated; T allele in rs 143383 was associated with 
the decrease of GFDS molecule expression and might 
increase susceptibility to osteoarthritis [17]. Moreover, it 
has been report that the differential binding of deformed 
epidermal autoregulatory factor 1 (DEAF-1) can modu- 
late the expression of GDFS via this polymorphism [29]. 
It is believed that GDFS plays role in the regulation of 
the chondrogenic cell growth and differentiation. These 
evidences support the function of GDFS gene which 
might indicate the importance of this polymorphism in 
osteoarthritis etiology. As the result, many scientists 
have been working for study the association of this poly- 
morphism in GDFS gene. Recently, a large scale associa- 
tion study between GDFS gene and osteoarthritis has 
revealed the genetic susceptibility of polymorphism in 

Table 3 Shows the association between allele and knee 



OA 



Allele 


Case 


Control 


Odds ratio (95%CI) 


p-value 


T 


117 


113 


1.53(1.013-2.306) 


0.04* 


C 


63 


93 


1 


N/A 



(C was used for calculation Odd ratio of T) 



Notes: ^statistically significant, N/A = not applicable 
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the GDFS; rsl43383 with odds ratio 1.12 (9596CI 0.99- 
1.31) [18]. This study included the data from diverse 
ethnics which have been reported in the recent litera- 
tures. However, the magnitude of this association is 
smaller than the Asian ethnic, it has been reported from 
Japan that odds ratio is 1.79 (P = 1.18 x 10-13) for 
GDFS polymorphism with per-risk allele (T) odds ratio 
by Miyamoto et al [17]. The magnitude of the associa- 
tion of GDFS polymorphism is different between Cauca- 
sian ethnic and Asian ethnic. Furthermore, the 
association of this polymorphism from Greek and Kor- 
ean study cannot demonstrate the association of GDFS 
polymorphism [20,30]. We have also reported the poly- 
morphism in ESR1, which is associated with knee OA in 
Korean population. However, our result did not show 
any statistically significant difference [31]. Thus, these 
evidences can imply that the association of ethic and 
genetic susceptibility in osteoarthritis might not be con- 
sistent among the different populations and cannot be 
applied to others. 

In addition, there is a report that GDFS polymorphism 
(rsl43383) also predisposes to Lumbar disc degeneration 
in women. Lumbar disc degeneration which is defined 
by disc space narrowing and the presence of osteophytes 
significantly associates with GDF5 polymorphism from 
cohorts from Northern Europe in women, with an odds 
ratio (OR) of 1.72 (95% CI 1.15-2.57) [32]. The genome- 
wide association study from Finland and Sardinia shows 
the common variants in GDF5 contribute to height dif- 
ference [33]. The GDF5 gene might be involved skeletal 
growth and development and GDF5 variants might play 
role in pathogenesis of bone and cartilage diseases. 

Although the sample size in our study was small when 
compared to previous studies, these limited samples 
were sufficient for demonstrating the statistically signifi- 
cant association of rsl43383 in GDFS core promotor 
area. Furthermore, the population in our study was in 
homogeneity and the distribution between the case and 
control groups was in the Hardy Weinberg's equili- 
brium. Therefore, our study can represent the significant 
susceptibility of GDFS in knee osteoarthritis in Thai eth- 
nic. Our findings have emphasized this association that 
the susceptibility of the polymorphism in GDFS among 
the Asian population because the magnitude of associa- 
tion is close to Japan and Chinese population. It is pos- 
sible that Thai population ethnic is believed to be close 
to Chinese ancestry. Finally, we decide to leave some 
suggestions in genetic susceptibility study. Firstly, it is 
important to study the genetic susceptibility in common 
disease as this might be valuable for further investiga- 
tion in order to understand the disease at molecular 
level or even apply to use for disease screening. More 
important, the association between polymorphism and 
disease might not similar in different ethics; therefore 



the genetic susceptibility from different ethics is 
required and the meta-analysis should be conducted in 
the adjacent area or similar genetic background. Lastly, 
other modern technology such as microarray technique 
should be employed for investigation of genetic suscept- 
ibility in osteoarthritis. 
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